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Abstract. We present the results of 5 BeppoSAX AOl 
Core Program observations of 3C 273 performed in Jan. 
1997 and compare them in detail with data obtained dur- 
ing the satellite Science Verification Phase (SVP), in Jul. 
1996 (Grandi et al. 1997). 

3C 273 was about 15% brighter in the first 1997 ob- 
servation than in the last one, and, on average, a factor 2 
brighter than the SVP observation. A count rate variation 
in the 2-10 keV band of ~ 12% in ^ 0.5 day was clearly 
detected during the last of the four pointings. 

Power-law fits with Galactic absorption to all obser- 
vations yield spectral indices in the range F = 1.53 — 1.6. 
Though a power law is an acceptable representation of the 
data in the whole 0.1-200 keV range, there is indication 
of a steepening of the spectrum as the energy increases. 
Residuals with respect to a single power law suggest the 
presence of a weak fluorescence iron line in the MECS 
data. No other features are detectable. Our data therefore 
mark a difference with respect to the SVP data, where a 
steeper power law below 0.5 keV, an absorption feature at 
~ 0.6 keV, and a more prominent fluorescence iron line 
have been found. 

The weakening of cold/warm matter signatures in our 
data with respect to the SVP ones may indicate that, at 
higher luminosities, the featureless continuum produced 
in a relativistic jet overwhelms any thermal and/or repro- 
cessed radiation, while the two components were almost 
comparable during the lower state of Jul. 1996. 

We quantitatively test this scenario, by considering an 
emission model which comprises the contribution from a 
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thermal Seyfert-like nucleus and a non-thermal power- 
law component, and find that indeed the observed fea- 
tures in both the AOl and SVP data are consistently 
reproduced by varying only the intensity of the non- 
thermal emission. Within this scenario, this radio-loud 
source shows evidence not only for thermal disk-like emis- 
sion, but also substantial reprocessing of X-rays onto cold 
matter. There is no evidence of a direct correlation be- 
tween the two components. 

Key words: Quasar: 3C 273; X-rays: galaxies 



1. Introduction 

3C 273 is a nearby (z = 0.158) quasar, and is one of the 
extragalactic objects best studied across the entire elec- 
tromagnetic spectrum. It shows almost all the features 
proper of high-luminosity quasars, i.e. radio and optical 
jet with high polarization, double radio lobes, superlumi- 
nal motion, variability at all frequencies, broad emission 
lines and signs of thermal emission in the UV band {Big 
Blue Bump). Multifrequency campaigns (e.g. Courvoisier 
et al. 1987) showed that the broad band Spectral En- 
ergy Distribution (SED) is rather complex, with two clear 
peaks at UV (^ 10 eV) and 7-ray (~ 1 - 10 MeV) ener- 
gies (Ramos et al. 1997, von Montigny et al. 1997). The 
general interpretation is that the radio-to-7— ray contin- 
uum is produced by the jet either via synchrotron and 
inverse Compton emission (e.g. Marscher & Travis 1996, 
Ghisellini et al. 1996), or by a proton-initiated cascade 
(Mannheim 1993), while the accretion disk is responsible 
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for the UV bump (see von Montigny et al. 1997, and dis- 
cussion therein) . Part of the accretion disk thermal emis- 
sion could be Comptonized in an accretion disk corona 
(e.g., Haardt & Maraschi 1991), or in the hot two temper- 
ature inner disk region, along the line originally proposed 
by Shapiro, Lightman & Eardly (1976), and more recently 
by Chakrabarty & Titarchuk (1995). 

The first detection in the X-ray band was obtained by 
Bowyer et al. (1970). Observations in the medium-hard 
X-ray range (up to 30 keV with GINGA) showed a hard 
power law continuum with photon index ranging between 
1.3 - 1.6 (Turner et al. 1990, 1991; Williams et al. 1992; 
Yaqoob et al. 1994), while both EXOSAT (Turner et al. 

1995) and, more recently, ROSAT data (Staubert 1992) 
reveal evidence of an excess above the extrapolation of 
the hard power law at energies ^ 1 keV. The study of 
the spectrum as seen by ROSAT (Leach, McHardy & Pa- 
padakis 1995) indicates the presence of a soft component 
(0.1-0.3 keV) with photon index F ~ 2.7 uncorrelated 
with the hard X-ray power law. Cold matter signatures 
such as a Compton reflection hump and an iron line were 
barely detected by GINGA (WiUiams et al. 1992). An 
emission line at 6.4 keV (source frame) seen at low flux 
levels was later confirmed by ASCA (Cappi & Matsuoka 

1996) . 

3C 273 shows prominent 7— ray emission as well, 
detected by OSSE (Johnson et al. 1995), COMPTEL 
(Hermsen et al. 1993) and EGRET (Hartman et al. 
1992), on board the Gompton Gamma-Ray Observatory 
[CGRO). OSSE data reveal that the hard power law ex- 
tends up to 1 MeV, with a break at higher energies 
(McNaron-Brown et al. 1995). 

Recently (Jul. 1996), the BeppoSAX SVP observation 
of 3C 273 has, for the first time, made available a simulta- 
neous coverage of the 0.2-200 keV energy band. Grandi et 
al. (1997) confirmed that the T ~ 1.5 power law extends 
at least up to 200 keV without any appreciable deviation, 
neither a Compton reflection hump nor a break. At ener- 
gies below 1 keV both an absorption feature and a soft ex- 
cess were present in the spectrum. A weak iron line at 6.4 
keV (source frame) was also present, with an equivalent 
width (EW) ~ 30 eV, consistent with GINGA (Williams 
et al. 1992) and ASCA (Cappi & Matsuoka 1996) resuhs. 

Here we present and discuss the results of pointings 
of 3C 273 performed in Jan. 1997. We will also make a 
comparison of these data with the SVP data, in view of 
understanding the long term behaviour of the source. The 
3C 273 observations discussed here are part of the AOl 
Core Program dedicated to bright blazars. 

The paper is organized as follows: in §2 observations 
and data reduction techniques are described, while §3 is 
devoted to the presentation of data analysis procedures. 
In §4 we present BATSE and optical data obtained dur- 
ing simultaneous observations. Results are then discussed, 
together with implications for theoretical models, in §5. 
Finally, a brief summary and conclusions are drawn in 



§6. Results of a preliminary analysis have been discussed 
in Maraschi, Fossati & Haardt (1997) and Haardt et al. 
(1997). If not otherwise indicated, hereinafter energies are 
measured in the observer frame. 

2. Observations and data reduction 

2.1. The BeppoSAX mission in brief 

For an exhaustive description of the Italian/Dutch 
BeppoSAX mission we refer to Boella et al. (1997). Here 
we briefly review the main characteristics of the instru- 
mentation, and the data reduction techniques applied. 

The narrow field coaligned instrumentation on 
BeppoSAX consists in a Low Energy Concentrator Spec- 
trometer (LECS), three Medium Energy Concentrator 
Spectrometers (MFCS), a High Pressure Gas Scintillation 
Proportional Counter (HPGSPC), and a Phoswich De- 
tector System (PDS). The LECS and MEGS have imag- 
ing capabilities in the 0.1-10 keV and 1.3-10 keV energy 
band, respectively, with energy resolution of ~ 8% at 6 
keV. At the same energy, the angular resolution is about 
1.2 arcmin (Half Power Radius). In the overlapping energy 
range the MEGS effective area (~ 150 cm^) is ~ 3 times 
that of the LECS. The HPGSPC covers the range 4-120 
keV, and the PDS the range 13-300 keV. In the overlap- 
ping energy interval, the PDS is more sensitive, while the 
HPGSPC has a better energy resolution. HPGSPC data 
will not be discussed in the present paper. 

2.2. The BeppoSAX observations 

The observations of 3C 273 were performed as part of the 
BeppoSAX AOl Core Program. The source was observed 
between Jan. 13-23, 1997. In this period 3C 273 has been 
observed 5 times: Jan. 13, observation period (OP) 1492, 
hereinafter referred as observation "A"; Jan. 15, OP 1502, 
observation "B"; Jan. 17-18, OP 1513, observation "C"; 
and Jan. 22-23, OP 1561 and OP 1563, observation "D". 
The datasets of the last two OPs have been merged into 
a single one due to the extremely short exposure time of 
OP 1561. The total effective exposure was 45.2 ks in the 
LECS, 92.1 ks in the MEGS, and 88.0 ks in the PDS. A 
journal of observations is given in Table I. For comparison, 
the net exposure times for LECS, MEGS and PDS in the 
SVP were 12, 131 and 128 ks, respectively. 

2.3. LEGS and MEGS data reduction 

Light curves and spectra for the LECS and the MEGS 
have been accumulated using the Xronos and Xselect 
tools (Ftools v4 . 0). We extracted the events from within 
a radius of 8' for the LECS and 6' for the composite MEGS 
image, consisting of the three co-added MEGS images. 
These selected events were then used to construct the 
light curves and to accumulate energy spectra for each 
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Table 1. Observation Log 



Obs. (start date) Net Exposure Time (ks) 
LEGS MEGS PDS^ 



Count rate (cts/s) 
LEGS'' MEGS'^ PDS'* 



A (01/13/97) 


13.6 


25.4 


22.8 


0.802 ± 0.008 


1.995 ± 0.009 


1.40 ± 0.04 


B (01/15/97) 


13.3 


24.1 


21.6 


0.729 ± 0.007 


1.883 ±0.009 


1.30 ±0.04 


G (01/17/97) 


12.5 


27.6 


25.0 


0.698 ± 0.008 


1.795 ±0.008 


1.24 ±0.04 


D (01/22/97) 


8.8 


22.5 


18.6 


0.695 ± 0.009 


1.691 ±0.009 


1.20 ± 0.04 



Exposure time for half of the PDS effective area. 
^ 0.12-4.0 keV count rate. 
" 1.6-10.5 keV count rate. 

13-200 keV count rate on half of the effective area. 



pointing. Different spectral binnings have been used and 
tested, such as the pre-constructed channel grouping tem- 
plates as provided by the BeppoSAX Science Data Center 
(SDC)Q or more standard groupings aimed to achieve a 
given minimum S/N across the entire spectrum. Unless 
otherwise indicated, the results of the spectral analysis will 
refer to spectra accumulated over the grouping2 template 
released by the SDC. 

LEGS data have been considered only in the range 
0.12-4 keV due to calibration problems at higher energies 
(Guainazzi 1997). The average count rate in the four 3G 
273 pointings is 0.726 cts/s in the LEGS, and 1.766 cts/s 
in the MEGS. Spectral analysis has been performed with 
the XSPEC 10 . GO package, using the response matrices re- 
leased in Sept. 1997. 

2.4- PDS data reduction 

The PDS was operated in collimator rocking mode, with 
a pair of units pointing at the source and the other pair 
pointing at the background, the two pairs switching on 
and off source every 96 seconds. 

Source visibility windows were selected following the 
criteria of no Earth occultation and high voltage stability 
during the exposures. In addition, the observations clos- 
est to the South Atlantic anomaly were discarded from the 
analysis. From collimators positions the ON and OFF time 
windows were also created and merged with the source vis- 
ibility window to create the final time windows on which 
the source+background and background spectra were ac- 
cumulated for each of the four PDS units, using the XAS 
software package, with an improved rise-time correction. 
The source was detected above 3 ct up to ~ 200 keV in 
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Fig. 1. LEGS hght curve. The bin width adopted is 12 ks. 



pointings A,B,G, and only up to ^ 150 keV in D. The 
grouped spectra from the four units were then coadded. 

In the calibrated energy band (~ 12-200 keV), the 
spectra were binned for all pointings in 15 energy inter- 
vals. The fits were performed with the XSPEC 10 . 00 pack- 
age using a response matrix appropriate for the improved, 
channel-dependent, subtraction of the background, for the 
correction for the dead layer of the PDS crystals, and for 
the improved rise-time correction. The average PDS count 
rate in the four pointings is 1.29 cts/s, on half of the ef- 
fective area. 



For 



details, 



see documentation 



at 



URLs 



http : / / www. sdc . asi . it /software /cookbook /spoctraL html _ 



and 



ftp://www.sdc.asi.it/pub/sax/cal/responsos/grouping/ 



3. Data analysis and results 
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Fig. 2. MECS light curve. The bin width adopted is 12 
ks. 

3.1. Temporal analysis 

In Fig. |l|, 1^ and || the overaU hght curves for the three 
instruments are shown. A decreasing trend is apparent in 
all the three instruments. We selected the MECS data for 
a detailed temporal analysis, as they are the BeppoSAX. 
data best suited for such a purpose, on account of better 
statistics and reliable performance stability. 

The MECS count rate monotonically decreases on time 
scale of days (Fig. The difference in the count rate 
between the first and last observation is ~ 15%. Similar 
variations are present also in the LECS and PDS data 
(Table 0), although the paucity of data does not allow 
to simply set limits on any lag between different spectral 
bands. We apphed the test to the total MECS light 
curve using time bin size ranging from 100 to 8000 s, and 
found that, independently of the time binning used, the 
count rate variation is significant at 99.99% level. 

The general trend shows up also in the LECS and PDS 
(see Table |, and Fig. |^ and ||), though the lower statistics 
does not allow a meaningful comparison with the MECS 
light curve. 

We also checked for MECS variability within single 
pointings. While the first three observations are consis- 
tent with a constant count rate, the last one is not, at 
the 99.88% significance level. The count rate reverses the 
decreasing trend and increases by about ~ 12% in about 
half a day. 

3.2. Spectral analysis 

We first present the broad band spectrum analysis and 
then consider each instrument separately. 
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Fig. 3. PDS hght curve. The bin width adopted is 12 ks. 

Table 2. Power Law** Fits: Broad Band 
Obs. T LE/ME'' PDS/ME'^ X^/d-o.f (xt) 



A 


+0.02 
^^•^'-0.02 


76+" "3 

'J- '"-0.03 


n SR+0 06 


63.2/84 (0.75) 


B 


1 cfi+0.02 


'J- '"-0.03 


Q4+0 07 


117.4/84 (1.41) 


C 


^•"^-0.02 


U.Oy_Q Q2 




92.4/83 (1.11) 


D 


J^-JO_0.03 


n 76+0 03 

"• '"-0.04 


0.831°:°? 


99.0/83 (1.19) 



Note: quoted errors are 90% confidence intervals for 3 

interesting parameters (Ax^ = 6.25). 

Nh fixed at Galactic value 1.69 x 10^° cm'^ 

LECS power law normalization relative to the MECS. 

^'^^ PDS power law normalization relative to the MECS. 

3.2.1. Broad Band Spectrum 

We jointly fit the data of the three instruments, for each 
of the four observations. The model adopted is a single 
power law with Galactic absorption (Nh = 1-69 x 10^" 
cm~^, Lockman & Savage 1995). We add two further pa- 
rameters, which we allow to be free, accounting for the 
mis-calibration of the LECS and PDS with respect to the 
MECS, which turned out to be the best calibrated instru- 
ment following independent tests. 

The data and the data/model ratio are plotted in Fig. |^ 
and Fig. ||, respectively, while the best fit parameters are 
listed in Table ^. As the fit to the broad band spectrum 
is dominated, in terms of statistics, by the MECS data 
points, the best fit value of the spectral index is essen- 
tially that of the MECS data (see next §3.2.3). As shown 
in Table ||, a single absorbed power law is an acceptable 
representation of the data for observations A and C, it is 
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top to bottom: A[xlOOO], B[xlOO], C[xlO], D 
1000 PT , — I — I ' — • — I 1 — I — 




0.1 1 10 100 

Energy [keV] 



Fig. 4. Data for the four observations. Observations A to 
D from top to bottom, respectively. For clarity, datasets 
A, B and C are scaled in count rate by 3, 2 and 1 orders 
of magnitude. 

marginally consistent with data set D, while it is rejected 
at 99.9% level in observation B. Visual inspection of Fig. ^ 
shows that large scatter is present in the PDS data, which 
also appear to be systematically steeper than the MEGS 
data. A feature around ~ 0.4 keV is present in the LEGS 
data. 

It is worth noting that, though the relative 
LEGS/MEGS and PDS/MEGS normahzations are 
marginally consistent with the fiducial ranges as given 
by the SDG (0.7-0.75 and 0.8-0.85, respectively), the 
LEGS/MEGS ratio possibly varies among different 
datasets. One should be careful in interpreting variations 
of the LEGS/MEGS normalization ratio, as the presence 
of the strongback and of a support grid structure in the 
LEGS window can actually affect such quantity (Parmar 
et al. 1997). Recently, an in-flight raster scan of the 
central (7' x 7') area of the LEGS and MEGS detectors 
(14 short observations of LMGXl) has allowed to map 
the LEGS/MEGS count ratio in the 2-8 keV band as a 
function of the position of the source in the LEGS field 
of view. The scan has shown that the LEGS/MEGS ratio 
can decrease by about 22% across the detector (Ricci 
1998, private comunication) . However, such effect occurs 
when the centroid of the source lies at least at 1.2' from 
the LEGS center. Since 3G 273 was always in the center 
of the field of view during the AOl observations, and the 
positions of the source centroid do not differ by more 
than 20", it is unlikely that window support obscuration 
affects the LEGS/MEGS ratio in our datasets. 

Nevertheless we must stress that fiducial values ap- 
ply only under the assumption that the fitted model is a 
good description of the data throughout all the three in- 
struments. If the model fails to describe the data in one 
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Fig. 5. Data/model ratio for the four observations. The 
fit model is a power law with Galactic absorption, with 
free normalization among the different instruments. 

of the instruments, this affects the resulting value of the 
inter-calibration parameter. 

We therefore conclude that a single power law, which 
is essentially determined by the MEGS statistics, is not a 
good description of the data in the LEGS and PDS instru- 
ment ranges. 

3.2.2. LEGS 

The LEGS spectra for the four observations were fit by 
a simple power law with free absorption. Data/model ra- 
tios are plotted in Fig. 6, while the best fit parameters are 
presented in Table ||. All the four data sets can be fit rea- 
sonably well by the absorbed power-law model. The 90% 
confidence regions for Nh and F are shown in Fig. |^. The 
spectral index is also plotted in Fig. |[ for Nh = Nn.gai- 
The value of F is in general smaller than that obtained 
in the broad band fit with fixed Galactic absorption. Ob- 
servations B and D require a value of Nh slightly lower 
than the Galactic absorption column density. This could 
indicate that a weak soft excess is present in the data, 
although the evidence is, at best, marginal. In fact the 
spectral indices for all the four observations are, within 
the errors, consistent with the average value F ~ 1.53 (see 
also Fig. 8). No apparent correlation of Nh with the flux 
is found. 

As evident from Fig. 5 and 6, in all the observations 
other than A, most of the contribution to the is due 
to residual structures at ^ 0.3 — 0.5 keV. It is known that 
they are spurious and produced by a mis-calibration of 
the LEGS in the Garbon edge region (Fossati & Haardt 
1997, Guainazzi & Grandi 1997, Orr et al. 1997). The 
latest LEGS response matrices coupled with the new data 
reduction pipeline strongly reduce the amplitude of these 
features, but do not completely remove them. Some excess 
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Table 3. Power Law Fits: LEGS 



Obs. 


Nh 

(10'° cm-2) 


r 


F[o.2-iikcv /d.o.f {xi) 
(10"^° erg/cm^/s) 


A 


1.69 (f) 

1.00_g 26 


J^-''*-0.03 
1 co+O.OS 


0.337 12.6/32 (0.40) 
0.336 11.3/31 (0.37) 


B 


1.69 (f) 

-I qr+0.27 
i.OO_Q 25 


1 c;7+0.02 
J^-0°-0.05 


0.341 42.2/32 (1.32) 
0.340 35.4/31 (1.14) 


C 


1.69 (f) 

'^■'"-0.30 


J^->J"-0.03 


0.332 25.1/32 (0.78) 
0.226 25.1/31 (0.81) 


D 


1.69 (f) 


1 t::4+0.03 

-, /(o+o.oe 

l-4»_o.o6 


0.219 37.4/32 (1.17) 
0.215 30.7/31 (0.99) 



Note: quoted errors are 90% confidence intervals for one (Ax = 2.706) or two interesting parameters (Ax = 4.61). 
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Fig. 6. Data/model ratio of power law fits to LEGS data. 



of residuals is also present near ^2.2 keV, probably due 
to inaccuracies in the calibration of the instrumental gold 
feature. 

Therefore, within these uncertainties, LEGS data are 
consistent with a power-law spectrum absorbed by Galac- 
tic material along the line of sight. In two out of four ob- 
servations (B and D) there is some marginal indication 
of the presence of a soft excess. There is no evidence of 
changes of the spectral index among the four observations. 

3.2.3. MEGS 

MEGS data of all of the four observations are well de- 
scribed by a single power law. Data/model ratios and best 
fit spectral indices are shown in Fig. 9 and 8 respectively, 
and their values are reported in Table ^. Gold gas absorp- 
tion column density was fixed at the appropriate LEGS 
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Fig. 7. 90% contour regions of column density and spec- 
tral index for absorbed power law fits to LEGS data. Solid, 
dashed, dot-dashed, and dotted lines represent observa- 
tions A, B, G, and D, respectively. The position of the 
labels within each region marks the best fit values. 



value in each data set. In any case, its precise value does 
not affect the fit of the MEGS data, as Nh is always lower 
than several times 10^° cm~^. 

With an average flux of ~ 1.14x10"^" erg/cm^/s [2-10 
keV] , the source was roughly in the middle of the range of 
fluxes historically observed (~ 0.6 — 1.7x10"^° erg/cm^/s) 
and the corresponding unabsorbed 2-10 keV luminosity in 
the Quasar frame is ^ 1.3 x lO^^ft^Q^ erg/s (Ji^q = i?o/50 
km/s/Mpc, go = 0). 

The energy spectral index is rather flat, with F ~ 1.55 
for three of the four observations (see Table consistent 
with the GINGA results (Lawson & Turner 1996), but 
slightly flatter than ASGA 1994 observations (Yaqoob et 
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Fig. 8. LEGS (filled squares), MEGS (open squares) and 
PDS (filled triangles) spectral indices for the four obser- 
vations. Error bars are 90% intervals for one interesting 
parameter. 



Table 4. Power Law'' Fits 


MEGS 




Obs. 


r F[, 

(10-1° 


-10] koV 

erg/cm^ /s) 


x'/d.o.f ixt) 


A 


1 ro+0.02 
1.00_Q 02 


1.18 


29.0/37 (0.78) 


B 


1 i;q+0.02 

1.00_Q Q2 


1.14 


42.1/37 (1.14) 


C 


1 64+002 
i.D4_Q 02 


1.07 


36.7/37 (0.99) 


D 




1.03 


48.7/37 (1.32) 



Note: quoted errors are 90% confidence intervals for one 

interesting parameter (Ax^ = 2.706). 

t'^' Nh fixed at LEGS best fit value (see Table |). 



al. 1994). The spectrum during observation C is instead 
steeper (F = 1.62±0.02). Given the small errors in the de- 
termination of the slope, the variation of spectral index, 
albeit small, is statistically significant. It is also impor- 
tant to note that the spectral index in the MEGS range is 
only barely consistent with the index derived for the lower 
LEGS energy range (see Fig. 8). 

The inclusion of an iron emission line (at 6.4 keV) does 
not statistically improve the fit in any of the four obser- 
vations. All the available MEGS counts could be added in 
a single spectrum to improve the statistics but, in view 
of the difference existing in the value of F between obser- 
vation G and the other three, we chose to sum the data 
normalized to the best fit of each observation, with the 
errors propagated accordingly. Furthermore, as weak nar- 
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Fig. 9. Data/model ratio of power law fits to MEGS data. 



Table 5 


. Power Law Fits: 


PDS 




Obs. 


F[20-1001keV 


x'/d.o.f ixt) 




(10-"' 


erg/cm'^ /s) 




A 


J^- "J-0.07 


2.58 


8.9/13 (0.69) 


B 


1 74+0.08 


2.62 


20.1/13 (1.55) 


C 


1 54+0.07 


2.32 


16.3/13 (1.25) 


D 


1 79+0.09 
^- ' ^-0.09 


2.17 


2.6/12 (0.22) 



Note: quoted errors are 90% confidence intervals for one 
interesting parameters (Ax^ = 2.706). 



row features could be smoothed out using a bin size that 
is too wide, in order to test for emission lines we binned 
the spectrum using the groupings template, which allows 
for a larger number of bins than grouping2. As shown in 
Fig. 10, both the shape of the residuals with respect to 
a constant (=1), and the sum of the ratios itself, clearly 
indicate the presence of an emission line at the source 
frame energy of 6.4 keV. The significance of such a fea- 
ture is above Scr, and its equivalent width is estimated to 
be ~ 20 eV. Although the residuals in Fig. |l^ might be 
even suggestive of the presence of an iron absorption edge 
at the expected energy, the statistics is not good enough 
to discriminate this feature from the noise. 

3.2.4. PDS 

Power law fits to the four PDS datasets alone are satisfac- 
tory in observations A, G and D (see Fig. |l^ and Table ||), 
despite a large scatter in the data points (probably related 
to a not-well calibrated background subtraction) . The ex- 
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Fig. 10. Data/model ratio (upper panel) and contribution 
to (lower panel) when the co-added MECS ratios (the 
data) are fit with a constant=l (the model). The dashed 
vertical line is at the energy of a 6.4 keV (source frame) 
iron emission line. 



ceedingly large value of the fit to observation B is es- 
sentially due to a single bin at ~ 100 keV, and, again, it 
is probably background related. The large error in the de- 
termination of the slope does not allow to detect spectral 
variations in the PDS data. The PDS spectral indices we 
derived are consistent with the OSSE results (McNaron- 
Brown et al. 1995). However, as shown in Fig. 8, they 
are statistically steeper than the MECS values in at least 
three out of four observations. 

We also analyzed the PDS spectra extracted using the 
SAXDAS software package, which, among other differences, 
does not includes the rise-time correction (see Section 
2.4). While the spectral indices are consistent with those 
reported in Table ^, the 90% confidence range is systemat- 
ically shifted to lower values for all the four observations, 
and marginally consistent with the MECS slopes. 

In conclusion, the steepening of the broad band spec- 
trum towards higher energies (Fig. ^) is evident using 
the pipelines of both of the two data reduction softwares, 
though the significance and degree of such a trend depends 
somewhat on the adopted PDS data reduction technique. 

3.3. SVP spectral analysis 

It is useful to compare the results of the BeppoSAX AOl 
observations with the previous BeppoSAX SVP data, pre- 
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Fig. 11. Data/model ratio of power law fit to PDS data. 



sented in Grandi et al. (1997) where the interested reader 
may find a detailed description of the observations and 
data reduction. 

We have repeated the LECS and MECS analysis of 
the SVP data taking advantage of the revised response 
matrices. The Grandi et al. results are, within the errors, 
quantitatively confirmed (see Table ^). 

The MECS continuum is described by a power law 
which is marginally steeper than that observed in Jan. 
1997, while the flux is a factor ~ 1.6 lower. MECS data 
also reveal a weak emission line interpreted as fluorescence 
Kq, transition of cold or mildly ionized iron. In order to 
estimate the line equivalent width, we operated as follows: 
first, we obtained the slope of the continuum by fitting the 
data excluding the 5-6 keV interval (model (e) in Table ^ , 
where the redshifted 6.4 keV line shows up; then we fit the 
continuum plus emission line keeping frozen the power law 
index, and the line energy and intrinsic width (model (f ) 
in Table |). 

In the LECS data a clear absorption feature and pos- 
sibly a soft excess below 0.5 keV, are superimposed to 
the power law. The fit to the 0.2-4 keV data with a bro- 
ken power law (the slope at energies ^ 0.55 keV is kept 
fixed at the MECS value) is not statistically distinguish- 
able from a fit with a single steep power law plus absorp- 
tion edge (see Table ^). However, as shown in Fig. |l^ the 
larger value of in the broken power law fit is entirely 
due to extra contributions occurring in the 0.5-0.7 keV 
range, i.e. there is a systematic defect at the energy of the 
absorption feature. 

To further check the reliability of the absorption fea- 
ture and soft excess detection, we computed the ratio be- 
tween the SVP and the AOl LECS data, as the ratio is 
certainly less subjected to instrumental effects. Fig. [l^ 
unambiguously indicates that the soft excess and the ab- 
sorption edge at ~ 0.6 keV are genuinely present in the 
SVP data only. 
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Table 6. SVP Data Spectral Fits 



LEGS 


Model 


Nh 


r 


break/ edge 

(keV) 


riiigh/ "^cdgc 


(10 


F[0.2-llkoV 

-1° erg/cmVs) 


x'/d.o.f ixi) 




PL'' 


1.69 (f) 


1 «9Q+0.036 
l.UZO_Q 








0.152 


43.2/32 (1.35) 




PT 


-0.31 


1 c;q7+0.068 

l.UO/ _Q Qgg 








0.148 


99 4/'Hl (0 95"! 




BKN PL'^ 


1.69 (f) 




61+°-^'' 


1.599 (f) 




0.154 


30.7/31 (0.99) 




PL+edge"^ 


1.69 (f) 


1 C71 +0.055 
i.ur l_o.060 


n 64+°" 

"•o^-o.io 






0.149 


24.0/30 (0.80) 


MEGS 


Model 


Nh 
(10^0 cm-2) 


r 


Ef. 

line 

(keV) 


(eV) 


(10 


F[2-10]keV 

-"-^ erg/cmVs) 


x7rf-o-,f (xi) 




PL° 


1.69 (f) 


1 500+0.010 








0.670 


41.1/32 (1.28) 




PL+line'' 


1.69 (f) 


1.599 (f) 


6.4 (f) 


on 9+5.7 
^»-^_11.6 




0.673 


46.9/37 (1.27) 



Note: quoted errors are 90% confidence intervals for one (Ax^ = 2.706), two (Ax'^ = 4.61), or three interesting parameters 



(Ax' = 6.25). 

single power law with Galactic absorption. 

single power law with free absorption. 
^'^^ broken power law with Galactic absorption. 

^'^^ single power law with Galactic absorption plus absorption edge. 

single power law with Galactic absorption, excluding the data in the range 5-6 keV. 
^'•^ single power law with Galactic absorption plus emission line. The line width is fixed at 0.1 keV, . 

Energy in the source frame. 



LECS, SVP: Broken Power Law 




LECS, SVP; Single Power Law + Edge 




0.1 0.2 0.5 1 2 

Energy [keV] 



Fig. 12. Contributions to ^ for fits to SVP LECS data: 
broken power law (top panel), and single power plus ab- 
sorption edge (bottom panel). 



Note also that the feature detected in the AOl 
LECS data at 0.4 keV vanishes almost completely in the 
SVP/AOl ratio. This implies that either it is equally 
present in both datasets or has an instrumental origin. 
We performed a survey in a sample of LECS datasets of 



different blazars (Fossati & Haardt 1997) and concluded 
that the latter option is the most likely. 

4. Simultaneous observations 

3C 273 has been constantly monitored by BATSE during 
the last few years, and it is of some interest to frame the 
BeppoSAX SVP and AOl observations into such a long 
look of the source. 

In Fig. H we show the BATSE light curve over 500 
days covering both of the BeppoSAX observations. The 
fluxes in the 20-100 keV energy band have been obtained 
by folding a single power law of photon index 1.7, after 
the PDS results. The BATSE data were collected by the 
Large Area Detectors (LADs) in Earth occultation mode. 
The occultation of a hard X-ray source in the sky by the 
Earth produces a step like feature superposed to the con- 
tinuous detector background count rate. A pair of rising 
and setting steps are generated, and the standard Earth 
occultation technique (Harmon et al. 1992) has been ap- 
plied to monitor X-ray sources on a daily basis. In the 
case of 3C 273 a bin of 20 days of the daily data has been 
applied in order to improve the statistics and enlight the 
trend of the source on longer time scales. 

Despite the fact that local minima or maxima are 
present in the BATSE light curve on monthly scale. Fig. [l^ 
shows that the BeppoSAX SVP observation was performed 
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comprises a 1242 x 1152 CCD camera with 0.47arcsec 
per pixel scale and standard Johnson-Cousins BVR fil- 
ters. Typical exposure times for 3C 273 are 180 s in the R 
band, 240 s in V, and 300 s in B. Data were analyzed with 
the Robin procedure locally developed, which performs 
bias subtraction, flat field correction and circular Gaussian 
fit after background subtraction. The source magnitude is 
obtained by relative photometry with respect to the ref- 
erence stars calibrated by Smith et al. (1985). Errors are 
of the order of a few hundredths of magnitude. Optical 
fluxes were calculated using the absolute flux calibration 
by Bessell (1979); no correction for galactic extinction is 
needed for 3C 273. 

The R, B and V average magnitudes are 12.7, 13.0 and 
13.1, respectively, corresponding to an integrated flux in 
R and V of - 1.2 x 10"!° erg/cmVs, of ~ 1.6 x 10"!" 
erg/cm^ /s in B. These values are slightly larger than the 
flux in the 2-10 keV band (see Table 4). It worths noting 
that no optical variability at a level larger than a few 
percent has been detected during the monitoring. 



Fig. 13. Ratio of SVP and AOl LECS data. Data/model 
ratios from observations A, B, C, and D of AOl are added 
together. 




10200 10400 10600 



TJD [days] 

Fig. 14. BATSE light curve. Twenty days intervals re- 
sulting from the rebinning of the daily data obtained with 
LADs in Earth Occultation Mode. The two vertical strips 
represent the windows of BeppoSKK observations during 
the SVP and the AOl. 



during a period of long term increase of the source, while 
the AOl observations occurred in a period of long term 
decrease. 

We also performed an optical monitoring during the 
AOl observations. The optical data were taken from Jan. 

12 to Jan. 31 at the Torino Astronomical Observatory, 
with the 1.05 m telescope REOSC. The instrumentation 



5. Discussion: The hidden Seyfert nucleus 

As reported in the previous section, the AOl data are 

only marginally consistent with a simple single power 
law model with Galactic absorption. There are (admit- 
tedly weak) indications of features, namely iron emission 
line and a high energy curvature of the BeppoSAX. spec- 
tra. More importantly, the presence of these same (but 
stronger) features, absorption around 0.6 keV, and a soft 
excess component at very soft energies diiring the weakest 
SVP state, lead us to consider a more complex underlying 
(variable) emission model. 

In particular, it has been already (qualitatively) sug- 
gested (see references in the Introduction, Cappi et al. 
1998) that an iron emission line in the spectra of radio- 
loud quasars might indicate that the X ray emission com- 
prises two components: a Seyfert-like thermal X-ray spec- 
trum, characterized by features induced by reflection over 
optically thick cold material subtending 2tt to the 'pri- 
mary' source, and a flat non-thermal component, arising 
as emission from plasma moving at relativistic speed in 
the jet, as typical of radio-loud objects. 

We therefore modelled the data with a complex model, 
consisting in the sum of a constant 'thermal' component 
and a variable non-thermal power law. Clearly our aim is 
not to obtain a significantly better fit for the AOl data, 
but rather to determine whether the variations in the ob- 
served features between the different observations are con- 
sistent with the scenario of an accreting object which dissi- 
pates energy in an accretion disc and its corona but where 
also a powerful jet is formed producing extremely variable 
beamed emission. 

The thermal spectrum has been represented by pa- 
rameters typical of Seyfert spectra (e.g. Nandra & Pounds 
1994) with a primary F = 1.9 power-law emission imping- 
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ing on cold matter subtending a 2tt solid angle, giving rise 
to an iron line at 6.4 keV (source frame) with an EW of 
150 eV. In order to establish the absolute normalization 
of the thermal component, we considered the SVP data. 
The flux of the underlying Seyfert-like spectrum is deter- 
mined by the simple requirement that the intrinsic 150 
eV of the iron line EW are observed, once diluted by the 
non-thermal component, as ~ 30 eV. This sets the flux of 
the Seyfert-like component at 6.4 keV (source frame) in 
the SVP to be 1/5 of the SVP total observed flux at this 
energy. We then assume that such flux level characterizes 
the Seyfert-like component in all the source states, i.e. the 
Seyfert-like component did not change among SVP and 
AOl data. In other words, we test the hypothesis that 
all of the spectral differences observed among SVP and 
AOl observations are accounted for by the variable non- 
thermal component. The latter has been modelled as a 
power-law with a (parabolic) turnover, to mimic the cut- 
off around a few MeV, detected in the 7-ray data (e.g. 
McNaron-Brown et al. 1995). 

Finally, as required by SVP-LECS data, we added to 
the model an absorption edge at low energy]^ making the 
assumption that it actually affects only the Seyfert-like 
continuum. 

We then fit this model to the LECS-MECS-PDS data, 
for SVP and the AOl brightest and faintest states (A 
and D). The three states are fit simultaneously, i.e. the 
edge energy and optical depth, the non-thermal power 
law index, and the LECS/MECS, PDS/MECS inter- 
calibrations constants are forced to be the same for the 
different sets of data in the fitting procedure. Only the 
normalization of the non-thermal component is allowed 
to be different for different datasets. 

The global fit is statistically good, with the best fit 
parameters reported in Table 0. The best fit model with 
all the different components relative to the SVP and AOl 
data is shown in Fig. |l^, and the data to model ratios for 
the three states are reported in Fig. |l^. 

The resulting values of the inter-calibration constants 
(LECS/MECS = 0.75 ± 0.01, PDS/MECS= 0.85 ± 0.03) 
are very close to the fiducial figures recommended by the 
BeppoSAX SDC. 

The edge is reasonably well constrained by the data at 
Eodgc = 0.62±0.08 keV with an optical depth t = l.Tlg J. 

The shape of the beamed component is determined by 
three parameters, namely the spectral index F, the energy 
above which the spectrum has a parabolic (in log) shape 
-E-break, and the vfi, peak of the spectrum E'pcak- 

The beamed component turns out to have F = 1.54 ± 
0.01, fully consistent with spectral indices detected in 
the brightest states of this source, when the non-thermal 
emission is likely to completely dominate the X-ray flux. 
In particular, we considered the recent results from ASCA 
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Fig. 15. Best fit model for the SVP and AOl-A data (top 
and bottom panel, respectively). The dashed and dash- 
dotted lines indicate the thermal reflected spectrum and 
the iron emission line, the dash-dotted one represents the 
non-thermal power-law emission and the solid line is their 
sum. 



by Cappi et al. (1998). The uncertainties in the calibration 
discussed by Cappi et al. do not allow a quantitative com- 
parison with the BeppoSAX results (see e.g. the discussion 
about the iron line EW by these authors). Nevertheless, 
we find that the model is also qualitatively consistent with 
the ASCA results, in particular with the hardening of the 
2-10 keV spectrum as the source gets brighter. 

The curved shape of the non-thermal continuum nicely 
accounts for the observed steep spectral index in the PDS 
range (see Fig. ||). 

The quality of the fit is rather insensitive to the energy 
of the peak of the parabolic component, and good fits 
are compatible with the power law starting steepening at 
energies above the MECS range and peaking around a few 
MeV. The values reported in Table are obtained with the 
peak of the parabolic component fixed at 2 MeV. 
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Table 7. Fit with Seyfert+Jet ModeP 
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Fig. 16. Data to model ratios for the SVP and AOl A and 

D datasets for the best fit parameters reported in Table ^. 



Parameter 



Value 



Units 



SVP 



#Nh 

LECS/MECS'' 
PDS/MECS'' 



General Parameters 

1.69 I 
0.158 
0.75 
0.85 



flO" 



+0.01 
0.01 
0.03 
0.03 



Seyfert Parameters 



#R 


1 




# Efoid 


100 


(keV) 


# cos(i) 


0.95 




#r 


1.9 




# Elinc 


6.4 


(keV) 


# Clinc 


0.1 


(keV) 


# Func 


3 X 10"^ 


(photons / cm^ /s) 


edge 




(keV) 


edge 


1 7+0.6 
-■-•'-0.7 


(10-1° erg/cm^/s) 


# F[2_l0keV] 


0.129 



Jet Parameters 



Our results show that BeppoSAX observations of 3C 
273 can be interpreted in a framework where the increase 
of the non-thermal emission in the A01~A swamps out 
the spectral features connected to the Seyfert-like com- 
ponent, such as the soft-excess, the absorption edge and 
the iron line. Moreover, it should be stressed that, in this 
model, the steep 'primary' thermal emission in the SVP 
state exceeds the non-thermal one at the lowest energies 
and naturally reproduces the steeper power law detected 
in the LEGS energy range. The steepening of the higher 
energy spectrum observed in the AOl data is simply inter- 
preted as the curvature of the inverse Compton emission 
reaching its peak. 

Once again, we stress that the model has not been 
constructed in order to best fit all the available data. The 
number of free parameters coupled with the data uncer- 
tainties would make this meaningless. Note in particular 
that there is no reason why any thermal component should 
be constant and there is actually clear evidence for vari- 
ability in the spectral slope (e.g., observation C in AOl), 
with no apparent correlation with the source intensity. 
Similar uncorrelated spectral variations have been also 
reported in EXOSAT and GINGA observations (Turner 
et al. 1990). Also, we can not exclude that the reflection 
hump is at a level lower than assumed (i.e., ri/27r < 1), as 
seems the case in Broad Emission Line Galaxies (Wozniak 
et al. 1998). 

However, from the fitting of the three states (SVP, 
low and high states of AOl) we conclude that Seyfert-like 
emission from an accretion disk and its corona is likely 
present, roughly at comparable flux level with the highly 
variable beamed non-thermal emission from a relativistic 
jet. The model in fact well reproduces three BeppoSAX 
observations and is consistent with the ASGA ones, by 



pb 

pb 
# Epeak 


1.54l°°} 
16.6l^y 
2.0 


(keV) 
(MeV) 


F[2-10keV] 

SVP; 


0.544 


(10-1° erg/cmVs) 


A: 


1.060 


(10-1° erg/cmVs) 
(10-1° erg/cmVs) 


D: 


0.892 


x'/d.o.f. ixi) 


230.6/238 (0.969) 



Note: 

Parameters marked with a have not been allowed to 
vary during the fit. The reported best-fit values and errors 
are those obtained keeping the value of the "Peak Energy" 
(Epeak) fixed at 2 MeV. The only parameter significantly 
affected by a change in Epcak is Ebreak- However for Epoak 
moving in the range 1 - 100 MeV it always holds 5 < 

Ebreak < 30 - 40 kcV. 

quoted errors are 90% confidence intervals for one 
interesting parameter (Ax^ = 2.706). 

this line flux gives an equivalent width EW = 29, 16 and 
19 eV for SVP, A and D datasets respectively. 



only allowing for a change in the normalization of the non- 
thermal emission. 



6. Summary and conclusions 

The results of BeppoSAX observations of 3G 273 per- 
formed in Jan. 1997 have been presented. The source var- 
ied within the observations and with respect to previous 
SVP data. Most notable is the variation of X-ray fea- 
tures clearly detected during the low state SVP, and barely 
present in the AOl data (and qualitatively consistent with 
the ASGA results recently reported by Gappi et al., 1998). 

The variable strength of the observed spectral features 
is consistent with a spectrum consisting of X-ray emis- 
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sion by a typical Scyfcrt- like nucleus, namely direct and 
reprocessed radiation dissipated in or above an optically 
thick accretion disk, and a variable non-thermal beamed 
component generated by the plasma moving in the rela- 
tivistic jet. Although the AOl data could be satisfacto- 
rily modelled by a single power law, slightly steepening 
above ^ 10 keV, the comparison with the SVP phase sug- 
gests that any thermal emission is indeed more evident in 
low source states (corresponding to weaker non-thermal 
radiation). Indeed the data available can be basically ac- 
counted for in this picture, by varying the absolute value 
of the normalization of the non-thermal emission. Clearly, 
monitoring of 3C 273 to derive variability information on 
the broad BeppoSAX energy range will allow to develop a 
more consistent and detailed model. 

Indications of thermal/Seyfert-like emission in 3C 273 
have been recognized for a long time, in different spectral 
bands, showing that in this source both thermal and non- 
thermal emission processes arc taking place, with com- 
parable radiative dissipation rate. The X-ray results fit 
well and support this picture. Further observations of the 
high energy broad band spectra of radio loud objects are 
critical in understanding the/any relationship between the 
dissipation in the accreting and ejected material. 
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